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A simple theory is developed that accounts for many of the observed physical properties of 
micelles, both globular and rod-like, and of bilayer vesicles composed of ionic or zwitterionic 
amphiphiles. The main point of departure from previous theories lies in the recognition and 
elucidation of the role of geometric constraints in self-assembly. The linking together of thermo- 
dynamics, interaction free energies and geometry results in a general framework which permits 
extension to more complicated self-assembly problems. 
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1. INTRODUCTION 

" Despite enormous progress in understanding the genetics and biochemistry of 
molecular synthesis we still have only primitive ideas of how linearly synthesised 
molecules form the multimolecular aggregates that are cellular structures. We 
assume that the physical forces acting between aggregates of molecules and between 
individual molecules should explain many of their associative properties; but 
available physical methods have been inadequate for measuring or computing these 
forces in solids or liquids." These few succinct opening sentences from the 
review of V. A. Parsegian 1 embrace and define a whole grey area bridging chemistry, 
physics and biology which is as yet but little explored. They imply a formidable 
injunction. For while it is axiomatic to the physicist or chemist that structural 
changes in any system should be reduced to a consideration of forces or free energies 
which cause those changes, the burden of proof lies with the proponent- The 
axioms of physics do not always receive so ready an acceptance from biologists whose 
whole thinking in the past has been centred on the role of geometry to the almost 
complete exclusion of forces and entropy. The burden of proof becomes especially 
great if one considers the increasing sophistication of those few successful theoretical 
advances in our understanding of condensed matter. To be convincing, and to have 
any hope whatever of reducing to some semblance of order the vast complexity of 
those intricate multimolecular structures that are the subject of biology, any successful 
theories of self-assembly must have as a minimal requirement extreme simplicity — 
to make them accessible to the biologist who has enough concerns of his own not to 
be dragged into the subtleties of modern physics. 

There is merit in the view that forces and entropy are important There is merit 
in the view that geometry is a determining factor in self-assembly. But there have 
been few attempts at model self-assembly problems which embrace both views. 
The aim of this paper is to develop a simple theory dealing with one area of self- 
assembly, the spontaneous aggregation of one-component lipid suspensions, where 
the fusion of both notions results in a unification of a multiplicity of diverse observa- 
tions. 

We shall deal with ionic and zwitterionic amphiphiles. These can take on a 
confusing variety of different shapes and sizes : some aggregate into small spherical 
or globular micelles, others appear to form long cylindrical micelles, while others 
coalesce spontaneously into vesicular or lamellar bilayers. 

The outline of the paper is as follows : In section 2 we delineate the necessary 
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thermodynamics which must underlie any description of self-assembly. Section 3 
is concerned with the formulation of a simple model, with minimal assumptions, for 
the free energies of amphiphiles in both aggregated and dispersed states. In section 4 
we carry out a preliminary application of the model to spherical micelles, and show 
that while the agreement between prediction and observation is qualitatively correct, 
the missing component can come only when geometric constraints are built into the 
model. Section 5 is devoted to a study of these geometric constraints. In section 6 
thermodynamic consequences of these constraints are examined in detail, and the 
refined model is shown to give remarkable agreement with observation. This 
confidence in our model having been established, we proceed in section 7 to study its 
consequences for the assembly of biological phospholipids which aggregate into 
bilayer structures. Finally we summarise our conclusions. 

2. THERMODYNAMICS OF SELF-ASSEMBLY 
(a) SOME GENERAL RESULTS 

To set notation and develop a framework for our later studies we briefly recapitu- 
late the thermodynamics of micelle formation. The literature is voluminous and 
confusing ; the most detailed statement of the problem being that of Hall and Pethica 2 
based on Hiirs small system thermodynamics. 3 We follow the approach of Tanford 4 
in which micelles (here meaning any lipid aggregate) made up of N amphiphiles are 
considered to be distinct species each characterised by its own self free energy JV/$. 
N will be called the aggregation number. Assume that dilute solution theory holds. 
Then equilibrium thermodynamics demands that the mole fraction X s of amphiphile 
incorporated into micelles of aggregation number N is given by 

kT 

li%+j r hi(X N IN)^ l i 0 l +kT\nX it (2.1) 

where k is Boltzmann's constant, T is temperature, and the suffix 1 denotes isolated 
amphiphiles. Alternatively 

X N = NX? exp mfi-Aytkll (2.2) 

For a prescribed total amphiphile concentration S, measured in mole fraction units, 
provided fi% are known functions of JVand T> the X N may be determined from eqn (2.2) 
together with the supplementary relation 

jg X n . (2.3) 

Two quantities which we shall require are the mean micelle aggregation number N 
and the standard deviation a of the distribution of sizes which measures polydispersity. 
These quantities are defined by 

<x J - <N-ff> 2 = £ N 2 X„f £ X N -N\ (2.5) 
It follows from eqn (2.2M2.4) that 

This equation relates the rate of change of micelle concentration with monomer 
concentration to the mean micelle aggregation number. Observed aggregation 
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numbers N tend to be large (typically > 40) since amphiphiles can decrease their free 
energies significantly by forming, and only by forming, large aggregates. Therefore, 
an immediate consequence of eqn (2.6) is that the micelle concentrations S~X t will 
be rapidly varying functions of monomer concentration X 1 . Define a critical micelle 
concentration (cmc) to be that value of X t for which S— X\ = Xf. For X x slightly 
less than the cmc X\ the micelle concentration will be much less than the monomer 
concentration, and for X x > X{ 9 the micelle concentration will be much greater than 
X t . Another useful expression which follows from eqn (2.3M2.5) is 

° ~B\nX x ~ dln(S-X x )' { 7) 

This equation relates the standard deviation to the rate of change of mean aggregation 
number with respect to micelle concentration (S- X x ). Clearly, a rapid change of 
with concentration is evidence of a large distribution of polydispersity in micelle size 
(at a given concentration). 

{b) CONSEQUENCES OF END EFFECTS 

The preceding relations are standard. If is negative, /.e., aggregates are 
energetically favoured, the transition from disaggregated to aggregated states comes 
about due to competition between entropy and energy, and little more can be said 
without spelling out the detailed form and magnitude of This will be 

reserved for later sections, but we anticipate here the forms which will emerge and 
make some general remarks. Mathematically only two possibilities exist for j$ 
which will allow the formation of large aggregates. We consider these in turn. 

(1) jtft may have a minimum value attained for some finite value of JV, say iV = M. 
This is the usual case discussed in the literature. 4 Physically, such a minimum or 
optimal value for p% comes about due to the competing effects, an increased hydro- 
phobic free energy of the hydrocarbon tails for N & M % and increased head group 
interaction due to electrostatic or geometric constraints for M. Very little can 
be said without a specific form of $ and a particular example of this case will be 
discussed extensively in section 4. A more convenient form for eqn (2.2) is here 

X fx \w u 

Two points can be made. First, even if $ = p% for all N > M, as will be shown 
below N cannot exceed M by very much for dilute lipid solutions. In fact N will be 
found in section 4 to be somewhat smaller than the optimal value M. Second, if the 
spread of micelle size is not large, the cmc will be given to a good approximation by 
neglecting polydispersity and assuming N = M. Whence (solving X x = X M for Xf) 9 
we have 

~G&-/*?)/*r. (2.9) 

This observation greatly facilitates analysis of experimental data. 

(2) The other possible form for $ is that p$ decreases with increasing # tending 
to a finite limit as N -» oo. This form will be seen later to arise from contributions 
to fx% from end effects, e.g., hemispherical-like ends on rod-shaped micelles, and 
hemicylindrical-like ends at the extremities of planar aggregates. 
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One functional form for $ which will occur in connection with rod-shaped 
micelles is 

= 00 l<N<m (2.10) 

where a is a constant, and m is a finite number below which micelles are energetically 
disfavoured, $ P & and, therefore, for convenience $ is set equal to infinity 
for 1 < N < nt. Substitution of this p% in eqn (2.2H2.3) leads to 



oo 



^ = e"« Y", S = X, +e"« £ NY", (2.11) 
where 

Y^X^xpivl-nDlkT. (2.12) 
The summation in eqn (2.11) can be carried out to give 

0 v mre-v y \ rtf .. 

s = ^+T3H 1+ ^Ty)> (2 ' 13) 

whence from eqn (2.5) 

For small 7, in fact for m(l-Y)P 1,N ~ m and as 7 increases to 1,# -> oo. But if 
m(l - 7) < 1 so that 7 « 1, we find 

(S-X t )e*~f^ 2 >m 2 . (2.15) 

Hence y»l corresponds to large and physically unreasonable micelle concentrations 
(S~X t ) unless a is large. That is, very large micelles can only occur at reasonable 
concentrations (typically in the range 10-M0- 2 mol dm- 3 ) if a is large. In this 
case 

This means that N is large and varies very rapidly with concentration S once the value 
of a exceeds about 25. In the other limit m(l - 7) > 1, we find 

(S-^e'-^. (2.17) 

This limit, therefore, corresponds to values of (S-X t ) e* small compared with m 2 . 
If we take m(l -Y) = x (with wj large and x ~ 1) we have 

_ m(s 2 +2s+2) 

~ x(x+l) 

(S-XJe'"^- (2.18) 



so that if x is chosen to make {S-X t ) e* ~ 1, we have N ~ ra. Therefore, with the 
form of /jft given by eqn (2.10), when a is small only small micelles (N ~ m) can occur 
at reasonable micelle concentrations (S— X x ) ^ 1. These micelles will be narrowly 
dispersed* When a is large, micelles will be large and broadly dispersed. 
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Another form which will be of interest in connection with the formation of 
bilayers is 

kT 

= oo 1 < N < m (119) 

which leads to the expression for the concentration 

S^Xi+f, NY N exp (-aiV*), (2.20) 

where 7 is given by eqn (2.12). This series has the property of converging on its 
circle of convergence (| 7| = 1), which leads to an interesting consequence. For a 
sufficiently large a, it will converge to a value of S < 1 for 7=1. Once this 
concentration is exceeded the amphiphiles will assemble spontaneously into infinite 
bilayers. (7=1 or $+kT\n X t = /& is equivalent to the condition that the 
chemical potential of amphiphiles in solution equals that of amphiphiles incorporated 
into infinite bilayers.) It is possible that the system will have a cmc before it reaches 
this phase transition. 

Consider as yet another possible example 

H° N = fx%+akT(N 2 ; N > m (2.21) 

which we shall use to illustrate the case when ix% vanishes faster than IjN. 
We find now 

S-X t = Y NY N e~*' N . (2.22) 

Since the result does not depend significantly on the lower limit of the sum when 
N > m we shall replace m by 1 for convenience. 
Consider the sum 

£ « £ Y» e-"" - i f *-m-im 9 (2.23) 

1 -oo 

where jS = -In 7. Applying the Poisson summation formula 5 to this, we find 

• 00 



N~~<X) J -oo 

= f W e-"*-"* dx+2 V f * cos {2nNx) e^""* dx. 

J 0 1 J 0 

By the substitution x = yj<x/p, this equation may be written 



(2.24) 



exp{-(oc^(y+l/^)}dy+ 

o 

CO 

2 Y I : z ! cos 



J 



[2*^)* yj exp {-Wfiy+m dy. (2.25) 



Consequently if £ = —In Yis small and a large, which is the limit of most interest 
to us, the period of the cosine in the integrand of eqn (2.2S) is much smaller than the 
width of the exponential term (the precise condition for this is <x*//J* > 1) and 
consequently the sum of eqn (2.25) may be neglected leaving 

E«d) i JJexp{-(a^)*(y+l/y)}dv. (2.26) 
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Let us further restrict ourselves to the case (a/?)* 1. The integral may then be 
evaluated by the method of steepest descents, 6 yielding 

Therefore 

(2.28) 

and 

8= —IniS-XJ-J^. (2.29) 
op 

This allows us to rewrite (2.28) as 

S-jr t ~^e- 2 «* (130) 

giving us an implicit expression for N as a function of S~X X . We have made the 
approximations (a#* ~ a/J7 > 1 and a*/£* = 27*/a > 1, £e., N <4*<NK One 
sees immediately from eqn (2.30) that under these conditions a large change in lipid 
concentration S—X x is required to produce a significant change in N, le. f N is a 
slowly varying function of 5- . In order that N may be a rapidly varying function 
of 5- X u it would be necessary that a ~ N, whereupon S- X x would be unphysically 
large. We conclude, therefore, that the form (2.21) would lead to lipid aggregates 
with little polydispersity. More generally, it can be shown that for a system with a 
chemical potential of the form p% = p% > +oJcT/N F > a phase transition to macroscopic 
aggregates will occur only if P < 1. For P > 1, the system will be narrowly 
dispersed about a finite aggregation number. When P « 1, the system will be 
broadly dispersed with a large aggregation number which is very sensitive to total 
amphiphile concentration. 

3. FORMULATION OF A MODEL FOR THE FREE ENERGIES 

No further progress can be made without some explicit model which gives the 
free energies p% 9 per amphiphile in both the aggregated and dispersed states. 
First consider the aggregated state. Contributions to the free energy p% fall into two 
classes, bulk and surface terms. It is well established that for most amphiphiles 
their hydrocarbon interiors in micelles and bilayers are in a liquid-like state above 
0 C C. 4 - 7 ' 1 1 Hence the bulk free energy per amphiphile, g % will be a function only of 
temperature T and the number of carbon atoms n. Surface contributions are of 
two kinds : (1) those arising from the attractive hydrophobic or surface tension forces, 
and (2) from opposing repulsive forces, in main of electrostatic origin. Since the 
hydrocarbon interior exists in a liquid-like state, we expect that the attractive force 
contribution can be well represented by an interfacial free energy per unit area of 
aggregate y, where y is close to 50 erg cm~ 2 characteristic of the liquid hydrocarbon- 
water interface. This value has been shown to be essentially the surface tension of 
water minus the dispersion energy contributions at the water-hydrocarbon inter- 
face. 12 - 13 As the surface tension of water varies by < 1 % over the range 0-0.5 mol 
dnr 3 NaCl, 14 we can expect the interfacial tension to be likewise insensitive to 
ionic strength. The surface area a per amphiphile is measured at the hydrocarbon- 
water interface, and the choice of this interface is a contentious issue. Herman 15 
and Tanford et al*> 16 choose this interface at a distance equal to the radius of one 
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water molecule beyond the van der Waals boundary of the outermost methylene 
groups. (Since dispersion contributions to y are negligible, this is equivalent to 
choosing the Gibbs dividing surface as the interface.) There appears to be no 
substantial reason for this choice. Indeed, thermodynamic arguments, 17 together 
with an analysis of some solubility data, 18 show that for various molecules the 
interface should be taken at, or even inside, the van der Waals boundary. We take 
the van der Waals boundary to define the interface. It will be seen subsequently 
that this choice goes some way towards removing a puzzling anomaly in our under- 
standing of hydrophobic interactions. We remark on one further problem in 
defining surface tension contributions ; The question as to whether we should take 
ya or y(a— a% where a' is the area per amphiphile covered by the head group and 
therefore not in direct contact with water. Since yd is constant, it can be absorbed 
into the bulk term g 9 to which it makes a small contribution. 

The repulsive surface terms are much more difficult to handle. The shape, size, 
and orientation of charged head groups, surface charge density, specific ionic 
adsorption, unknown dielectric constant of the surface region, certain occurrence of 
Stern layers and associated discreteness of charge effects all conspire to inhibit any 
rigorous analysis. For example, in recent work on ionic micelles, Stigter 11 showed 
that use of a Chapman-Gouy model for the double layer based on the non-linear 
Poisson-Boltzmann equation gives an outward electrostatic pressure which varies 
by only 20 % over an electrolyte concentration of 0-0.4 mol dm- 3 NaCl. On the 
other hand, the cruder Debye-Hiickel approximation taken with an arbitrary factor 
of the order of \ gives a reasonable description of the critical micelle concentration and 
other properties ofionic lipid micelles. 4 * 19 » 20 As regards the Stem layer, Stigter 11 • 21 
has carried out an analysis for a variety of ionic micelles and concludes that the 
counterions are effectively in a layer of thickness 4-5 A. Despite the apparent 
intractability of the problem, the various approaches suggest that all of these 
complications can be subsumed by simple phenomenological forms. Various 
authors since Debye 22 have attempted such an analysis. 11 ' 19 " 26 Thus, a repulsive 
energy contribution which varies as a constant/a has been shown by Tanford 19 to 
give a realistic description of micelle size and cmc. A repulsive energy of this form 
would arise from a double layer of charge as in a capacitor, with charge e(a per unit 
area, separation D of the capacitor planes, and dielectric constant e. The magnitude 
of the constant is then 2ite 2 Dfe. We adopt this form. 

We stress that this choice for the repulsive energy is not critical. All the sub- 
sequent analysis can be carried out with other repulsive forms, e.g., constant/a 2 . 
The advantage of the capacitance model is that it gives about the simplest expression 
for the energy with any semblance of physical reality which is amenable to exact 
analysis. Such an assumption clearly disguises a multitude of sins, and can only be 
justified a posteriori. 

One further advantage of this form is that for zwitterionic amphiphiles where 
the head groups are normal to the surface (e.g., as is believed to be the case for 
lecithin 21 * 28 ), the value of D can be readily identified with zwitterionic charge 
group separation. The expression we have chosen is to be expected here (see 
Appendix A). Additional support for the. universal phenomenological form constant/ 
a is furnished by the observation 29 that both ionic and zwitterionic amphiphiles 
have very similar micellar properties. If these assumptions are granted, the free 
energy /# per amphiphile in the aggregate is 

fl ° N =ya + +g. (3.1) 
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(a) CURVATURE CORRECTIONS 

The form eqn (3.1) is to hold strictly for planar surfaces. We have already 
alluded to curvature corrections to the effective surface tension contribution. If the 
micellar surface is spherical, the effective surface tension y for both a drop and a hole 
in the liquid is less than that for a planar interface by a factor (1-8/R), where <5 
is of the order of one or two molecular radii and R is the position of the Gibbs 
dividing surface. 17, 18 Corresponding corrections to the electrostatic contributions 
are expected to be of much more importance, and can be handled within the frame- 
work of a capacitance description. Thus for a spherical capacitance the energy per 
unit area is from electrostatics 

2ne 2 D 

£ ' " eaO+DIR)' ° 2) 

Except for two instances in the subsequent development, curvature corrections can 
be ignored. 

(b) MONOMER FREE ENERGY 

The monomer free energy comprises two main contributions, a hydrophobic 
term g' due to the self-energy of the hydrocarbon chain in water, and an electrostatic 
self-energy associated with the (charged) head group. The term g' is undefined by 
itself, but the important quantity g'~ g has a meaning as the usual hydrophobic 
energy 4 required to take a hydrocarbon chain from water to bulk hydrocarbon. 
For this quantity we use the measured value so that a precise model for this contribu- 
tion is unnecessary. Characteristically * g'-g ~ 825 cal mol -1 (CH 2 group)- 1 and 
2100 cal mol- 1 (CH 3 group)- 1 . (We note in passing that in analysing the measured 
hydrophobic free energy changes for n-alkanes, various authors x5, 16 assign a value 
to the hydrophobic interfacial energy of y h = 20-33 cal moH A~ 2 , equivalent to 
14-23 erg cm -2 . The range of values reflects the particular choice of " cratic " 
contributions. 4 The assignment of such values depends critically on choice 
of the surface of tension, normally taken at the Gibbs dividing surface. As 
already remarked, this is probably incorrect, and if the vander Waals boundary of 
the hydrocarbon is chosen, a higher value close to the bulk hydrocarbon-water 
interfacial energy obtains, i.e., 50 erg cm -2 ,) 

Similarly, the change in electrostatic self free energy of the head group between 
aggregated and dispersed state is the relevant entity required. As shown in Appendix 
A, the capacitance form for the repulsive surface term includes only head group- 
head group, head group-counterion and counterion-counterion interactions ; Le. 9 we 
have chosen a reference state in which the head group self-energy is omitted. 
Although the head group in the aggregated state may have a different self-energy 
from that in the dispersed state due to its proximity to a region of low dielectric 
constant, both free energies, and therefore their difference, can be expected to be 
much smaller than the hydrophobic free energy change. (A very crude estimate of 
the electrostatic self-energy is provided by the Born expression e 2 /2er. For e « 80, 
r«2A, we have e 2 /2er « 1.7 kT while typically hydrophobic contributions to 
/# amount to £ 15 AT.) Hence we ignore changes in electrostatic self-energies. 
The interaction free energy of a head group in the dispersed state with counterions 
and other head groups is certainly negligible, as is evident fiom measured activity 
coefficients of salt solutions. 30 

For zwitterionic amphiphiles there will be a small additional electrostatic 
contribution to of the form £ 2 /(e|r + -r-|). 
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This completes the formulation of our model for micellar free energies, except for 
one caveat : to examine now the consequences of the model, the thermodynamic 
formulation of section 2 requires that the relation between free energy (a function of a) 
and the aggregation number N be spelt out. The connection between a and N is 
determined by the shape and size of the aggregate. For a given a y fixed hydrocarbon 
volume per amphiphile and aggregation number, there are geometric or packing 
constraints which specify the allowed shape and size. An examination of these 
packing constraints and their consequences will be deferred to sections 5 and 6. 
For the moment we shall be concerned with testing the model for the simpler problem 
of spherical micelles, and comparing with experiment, in order to develop confidence 
that the assumptions embodied in eqn (3.1) provide a satisfactory basis for prediction. 

4. PRELIMINARY APPLICATION TO SPHERICAL MICELLES 

To illustrate and test the viability of the model, we now consider its application 
to the formation of spherical micelles. In this application we ignore the important 
geometric constraints dealt with in the following section. These complications will 
be built into the theory in section 6. We consider spheres for the following reasons. 
If we ignore packing constraints, of all the possible shapes for aggregates, the sphere 
is the smallest shape with a given surface area per molecule and therefore has the 
lowest aggregation number. Consequently a sphere will be the thermodynamically 
favoured state, even in the absence of curvature corrections to electrostatic free 
energies, or of end-effects which will occur for rod-like micelles or planar bilayers. 
If curvature and end-effects are included, non-spherical shapes, e.g., cylinders and 
bilayers, turn out to be even less favoured. Although spherical micelles are unlikely 
to occur in general due to packing constraints, 29 * 31 restriction of the discussion to 
these aggregates gives us an idea of the behaviour to be expected of globular aggregates 
as opposed to large cylinders or bilayers. Distortions from spherical shape change 
the relation between area per amphiphile a and aggregation number N and this 
modification can be incorporated subsequently. 

The thermodynamic equations of section 2 are expressed in terms of N 9 whereas 
given in section 3 involves a. We need a relation between a and N. The 
geometry of a spherical micelle is illustrated in fig. 1 and if v denotes the volume of 
the hydrocarbon tail of an amphiphile we have the relations 

and 

The first relation assumes that the density of hydrocarbon is constant in the interior 
of the micelle. 11 If we ignore curvature corrections the free energy per amphiphile 
in a micelle is from eqn (3.1) and (4.2) 



A 



= \*+ 7 )+il Co - . (4-3) 

Henceforth a 0 will be referred to as the optimal surface area per amphiphile, being 
that area at which the free energy per amphiphile in a micelle is a minimum. From 
eqn (2.8) the distribution of aggregates is given by 

X„ = Jvjj£ exp W (4.4) 
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Take M to be the aggregation number at which p% takes its optimal value, /.&, 

&**2ya 0 +g, M«^p-, = T (4.5) 

so that eqn (4.4) takes the form 

Note that since a 0 = (47i)*(3i;)*/Af * and the tail volume of an amphiphile is known, 
the only parameters which occur in the distribution eqn (4.6) are X M% M and y, the 
interfacial tension of the hydrocarbon-water interface. 



o 




• o 



• 

Fig. 1. — Spherical micelle. 




N 

Fig. Z— Concentration of amphiphiles X* in spherical micelles of aggregation number N. The 
curves are plotted for three assumed values of dielectric constant e. The curve for e = 0 is equivalent 
to ignoring curvature corrections to the electrostatic energy. The curves are normalised to the same 
maximum value. Xm is taken as 10~ 4 mole fraction (5 mmol dnr 3 ). 

Fig. 2 shows a plot of X H against N for M = 60 and M = 100. For these plots v 
was taken to be 350 A 3 corresponding to a hydrocarbon chain of 1 2 carbon atoms, 4, 1 1 
and y = 50 erg cm* 2 . Several features deserve comment. The maximum of the 
distribution is insensitive to X M , and therefore to total lipid concentration. The 
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curves are essentially Gaussians distributed about the maximum so that the mean 
aggregation number N and the maximum are identical. N is reduced from the 
optimal value M by about 10 %, so that in view of eqn (4.2) and (4.5), the area a at 
the maximum is only about 4 % greater than a 0 , or less for higher values of Af. 
The standard deviations of the curves are given by 

and the spread about the maximum is therefore a narrow one. Also, as expected 
from eqn (2.7) a change in lipid concentration X u by an order of magnitude changes 
the mean aggregation number by only 1 or 2. These distributions are similar to 
those derived by Tanford. 20 The distribution does not depend on the specific form 
for Any admissible form of fi N which comes about due to a balance between 
opposing forces would lead to a similar distribution. 

In comparing with experiment, the actual distributions are not known and are 
difficult to determine. What is known and available is the cmc under a range of 
conditions. In general, precise determinations of In (cmc) for a given choice of 
parameters involve a complicated numerical routine. However, as already remarked, 
the system is narrowly dispersed, and can be treated to a good approximation as a 
monodisperse system. Further, the area a corresponding to the mean aggregation 
number N differs by only 4 % from the optimal area a 0 . If curvature corrections to 
H% are included a is shifted even closer to a 0 . Thus eqn (4.3) becomes 

-<° + ,-pfy + '- < 4 - 8) 

Since DfR = 6vEyfMe\MfN)*, after a little algebra, repetition of our previous 
analysis gives, instead of eqn (4.6), a more complicated algebraic expression involving 
an additional parameter e. The corresponding distribution is plotted in fig. 2 for 
e — 80 and 40, and we see that the net effect of curvature corrections is to shift N 
towards M and a towards a 0 - Thus, with little error, we can take S = a 0 > R - 
and use eqn (2.9) which becomes 

In (cmc) ~ (2ya 0 +g-g')lkT 




(4.9) 



Before we put numbers into these equations we note some qualitative predictions. 
Since v is roughly proportional to the number n of carbon atoms in the hydrocarbon 
tail, eqn (4.5) shows that (K ~ JR M ) is proportional to n as observed. 32 Further, 
q 0 is independent of n. This is consistent with experiments on monolayers at the 
oil-water interface where the pressure-area isotherms are found to be insensitive to 
chain length above the transition area. 33 As g 1 — g is a linear function of n, In (cmc) 
is also linear in n 9 again as observed. 4 

The much studied sodium dodecyl sulphate (SDS) micelle provides a more 
substantial test. Huisman 34 measured the variation of cmc and N for these micelles 
as a function of NaCl concentration from 0 to 0.3 mol dncr 3 (see also Emerson and 
Holtzer 35 ). In eqn (4.9) all parameters except M are weakly dependent on salt 
concentration. Fig. 3 shows the plot of In (cmc) against M~*. The best fit straight 
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line gives K x = 44, K 2 = 20. With r S i>s = 350 A 3 , 4 K x » 44 implies that y ~ 37 
erg cm- 2 . This value, while reasonable, is somewhat less than the value 50 dyn cm*" 2 
for the oil-water interface. This empirical value for y allows an estimate to be made 
of the capacitance distance Z>. Taking a = 80, D varies in the range 8 to 5 A as salt 
concentration varies from 0 to 0.3 mol dm*" 3 . For smaller values of i the values of 
D would be proportionately smaller. We do not place much emphasis on the magni- 
tude of this parameter, but note that it is of the order of magnitude one would expect 
for the size of the Stern layer. 11 Next consider the value of (g'—g) = K 2 kT = 
20&r= 12 000 cal mol*- 1 . This compares with the expected hydrophobic free 
energy for transport of 1 CH 3 + 11 CH 2 groups from water to bulk hydrocarbon of 
~11 000-12000 cal mol- 1 . 4 If we make allowance for the area d covered by the 
head group, the best fit hydrophobic energy would be reduced to ~ 10 500 cal mol" 1 . 
If in addition we admit curvature corrections this value is reduced further to ~ 10 000 
cal mol" 1 . 



in 




20 .21 .22 .23 2L .25 .26 .27 



Flo. 3.— Variation of In (cmc/55.5) with M~* for SDS in NaQ solutions as given by eqn (4.9). 
Points : experimental results. Theoretical curve is based on the assumption that the micelles are 
spherical, cmc/55 J is in mole fraction units. 



We could go further at this stage and check the predicted n dependence, effect of 
temperature on cmc, and analyse changes with salt in more detail All of these are 
in good semi-quantitative agreement with experiment To sum up : the model as 
developed so far gives good agreement with observation, although it might be con- 
sidered that the surface free energy appears 30 % too low. However, as we cannot 
yet explain the existence of rod-like micelles or bilayers which suggests that geometric 
constraints cannot be ignored. We postpone a more precise analysis until the effects 
of these constraints have been investigated. 



5. PACKING CONSTRAINTS AND NON-SPHERICAL MICELLES 

We have already alluded to geometric limitations which place restrictions on the 
allowed shapes of micelles, and it is clear that packing constraints must be invoked 
for a proper treatment of self-assembiy, for, in the absence of any such restrictions, 
spherical micelles will always be thermodynamically favoured over other shapes like 
cylindrical micelles or bilayers. There must then be some overriding factor that 

IT— ^9 
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forces some amphiphiles to assemble into these larger structures which appear to be 
thennodynamically disfavoured. 

In this section therefore, we study the packing properties of amphiphiles which 
come about due to geometry. The most general form of the problem is exceptionally 
tedious because thermodynamics and packing are inextricably linked and cannot be 
considered in isolation. Further, the general problem of relating aggregation 
number to shape involves a difficult problem in the calculus of variations. Nonethe- 
less we can make considerable progress in elucidating the role of packing provided we 
assume that the mean area per amphiphile wilt always be close to the optimum area 
a 0 . As shown in section 4, this assumption is fair for systems which exhibit little 
polydispersity, and will allow us to disentangle energetic and geometric factors from 
the thermodynamics. This leads us to consider the shape and size of the smallest 
micelle consistent with a = a 0 . This we call the critical packing shape. We shall 
find that beyond some characteristic critical aggregation number N c , both spherical 
and ellipsoidal micelles are prohibited. Geometric factors lead us to consider 
different shapes once these are excluded. 

Consider first a spherical micelle (fig* 1). Here the radius R, hydrocarbon core 
volume v and surface area a per amphiphile at the hydrocarbon-water interface are 
related by 

3v 

- = (5.1) 
a 

Then, since the radius of a spherical micelle cannot exceed a certain critical length, l c , 
roughly equal to but less than the fully extended length of the hydrocarbon chain, it is 
apparent from eqn (5.1) that once vfa 0 l c > J, spherical micelles will not be able to 
form unless the surface area a > a 0 . This gives a critical condition for the formation 
of spheres : 

On the assumption that the surface area per amphiphile is everywhere equal to or 
close to the optimum area a 0 , we must therefore look for alternative non-spherical 
shapes once v(a 0 l Q > |. (The way v and / c are actually determined from the nature 
of the hydrocarbon chains is model-dependent and not of immediate concern here.) 
Similarly, the critical condition for cylindrical micelles and planar bilayers is readily 
found to be respectively 

v 1 v 
— — -L (5.3) 
a 0 l c 2 a 0 l € 

The transition shapes of amphiphilic structures, as they go from spheres to 
cylinders 36t 37 have not been studied ; this will be our next concern. 

Regardless of shape, any aggregated structure must satisfy the following two 
criteria : (1) no point within the structure can be farther from the hydrocarbon- 
water surface than / c ; (2) the total hydrocarbon core volume of the structure V and 
the total surface area A must satisfy V/v = A/a 0 = N. This criterion is only 
approximate since it assumes that the mean surface area per amphiphile is equal to Qq. 

(a) ELLIPSOIDAL MICELLES 

Tarter 31 and more recently Tanford 29 have examined the way amphiphiles can 
pack into prolate and oblate spheroids once they can no longer pack into spheres. 
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The shorter ellipsoid semi-axis was put equal to / c and the longer axis was then 
determined by applying criterion (2). Both oblate and prolate ellipsoids are found 
to satisfy criterion (2), though for a given v, l 6 and a 0 oblate spheroids have a slightly 
lower aggregation number N. Clearly, criteria (1) and (2) can be satisfied by many 
different shapes. Closer scrutiny of spheroids shows that these shapes are in fact 
untenable and leads us to a stronger packing condition. This can be seen as follows. 




W (6) 
Fig. 4.— Oblate spheroids. 

Consider an oblate spheroid of large eccentricity [fig. 4(a)] satisfying criterion (2) 
and with minor semi-axis equal to / c . Any local region of the oblate spheroid looks 
very much like a bilayer. However, for the amphiphiles to be so packed we require 
that at = t? (a constant), where It is the (varying) hydrocarbon thickness. Conse- 
quently the area per amphiphile satisfies a ^ a 0 almost everywhere so that the 
amphiphiles of the spheroid must also be energetically unfavourable almost every- 
where. [Any form for the opposing surface forces will give the surface energy per 
unit area as constant xa 0 +consta.ntx(a-a 0 ) 2 .] Criterion (2) disguises this 
energetically unfavourable situation since it contains the volume V and area A of the 
whole spheroid, found by averaging the amphiphile volumes and surface areas over 
the whole structure. By averaging the volumes and areas in this way the discrepancy 
in the local surface area and energy per amphiphile is averaged to zero, even though 
the local packing is almost everywhere energetically unfavourable. A similar 
situation holds for a prolate spheroid of large eccentricity. Thus criteria (1) and (2) 
are necessary but insufficient conditions for viable micellar structures. Another 
criterion is required which would indicate whether amphiphiles are packing in an 
energetically favourable way at all parts of an aggregate and not only on the average. 

(b) LOCAL PACKING CRITERION 

Consider an amphiphile (fig. 5) of surface area a and liquid hydrocarbon core 
volume v in a micelle (or bilayer vesicle) where the local radii of curvature are R x and 
/? 2 . Referring to fig. 5 we obtain by elementary geometry a " packing equation 99 

where / is the length of the hydrocarbon region of the amphiphile. This equation is 
exact for spherical surfaces (J? x = R 2 \ cylindrical surfaces (R 2 = oo) and planar 
surfaces (R t « R 2 « oo), and holds to a high degree of approximation, with an 
error of not more than 1 %, for surfaces of arbitrary curvature. 

The weaker criteria (1) and (2) are embodied in eqn (5.4) if we stipulate that 
/ < / c . Application of eqn (5.4) to any part of a micelle or bilayer provides a test for 
whether or not amphiphiles there can pack in an energetically favourable way, i.e., 
with a = a D . There are no 14 new concepts " here ; the packing equation is simply a 
necessary geometric expression that relates t>, a y /, R t and 
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To illustrate the application of the packing equation we now show that spheroids 
cannot accommodate amphiphiles in their regions of greatest curvature once these 
amphiphiles can no longer pack into spheres, ue. 9 once &/064 exceeds Consider 
the oblate spheroid of fig. 4(A). Since by assumption spheres are prohibited, ^ 
have r t > / c , and by criterion (2), r 2 < / c . In the region of greatest curvature />, 




Fig. 5. — Geometric packing of a hydrocarbon region of volume v and surface area a at a surface with 
two radii of curvature Ri and R 2 - Eqn (5.4) gives the relation between v, a, R u Ri and the length 
of the hydrocarbon /. For both R u R 2 > /, a void region is formed behind the hydrocarbon region. 

the principal radii of curvature are JR t = r u and R 2 = r 2 /r lf so that R 2 < / c . Since 
vfa 0 > / c /3, it is sufficient to show that v[a 0 as given by the packing equation is less 
than / c /3. We have from eqn (5.4) 

since / must satisfy / < R 2 . Successive application of the inequalities J? 4 > R 2 > 
R 2 < l c gives further v/a < R 2 /3 < / c /3 which establishes the result. Thus the 
oblate spheroid, rather than being a viable alternative shape when amphiphiles can 
no longer pack into spheres, is seen to be energetically unfavourable the moment the 
spherical shape itself becomes unfavourable. Indeed, even when amphiphiles can 
pack comfortably into spheres, eqn (5.4) shows that only spheroids of low eccentricities 
can accommodate the amphiphiles as well. Likewise for prolate spheroids. 

(6) PROPERTIES OP ALLOWED SHAPES 

Strictly, the packing equation should be all that is needed for a complete solution 
of the problem of finding the shape (or shapes) that can pack amphiphiles of a given v t 
a 0 and / c . However, as already remarked, the exact solution to the problem is both 
difficult and uncalled for. On reflection, it is apparent that the simple shape which 
satisfies all the packing conditions is not unlike a distorted oblate spheroid. Recall 
that the unacceptability of an oblate spheroid comes about because the peripheral 
regions had too great a curvature while the central regions were too thick. If> 
however, the central region is flattened so that locally it approaches a bilayer and 
simultaneously the curvature of the peripheral regions is reduced, the packing criteria 
can be satisfied. Of a number of shapes that have the general features of such a 
distorted oblate spheroid, a simple analytically tractable shape which satisfies all 
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our criteria is illustrated in fig. 7. This shape may be termed an ellipse of revolution. 
Qualitatively, it has the following properties : when amphiphiles can still pack into 
spheres the ellipse is a circle centred at the origin and the shape is therefore a sphere ; 
as v/a 0 increases above IJ3 the circle distorts into an ellipse of small eccentricity 




Fig. 6.— Wedge-shaped volume take up by hydrocarbon when R 2 = L Note that the void region 
existing when 2?i > / and R 2 > I has disappeared (cf. fig. 5). 



whose centre moves progressively outwards from the symmetry axis, and the micelle 
becomes globular. If this process were continued indefinitely, at some value of 
vjciQ the shape would become toroidal-like. Ultimately, as v/a 0 approached / c /2, 
the critical condition for cylinders, the ellipse would have become a circle once again 
giving an infinite cylinder. 

sphere globule toroid infinite cylinder 
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Fig. 7.— Approximate transition shapes from spheres to cylinders, based on purely geometric con- 
siderations. For any given shape the area per amphiphfle a 0 is uniform. The minimum aggregation 
number M is achieved at / = / c as explained in the text. More realistic transition shapes may be 
different due to thermodynamic requirements, discussed in section 6. 

In the globular regime, ue., when micelles are still near spherical, we expect, as 
shown in section 4, that the distribution of sizes will be narrow, so that the approxi- 
mation a ~ a 0 is good ; any shape which satisfies the packing criteria will be adequate 
here. When we come to larger aggregation numbers, we are on much more tentative 
ground, because the assumption of monodispersity is no longer valid, and the average 
micelles can be far from optimal. Indeed, when thermodynamic considerations are 
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involved, the initial growth phase discussed in section 6 turns out to be by a process 
which is in fact energetically unfavourable, through rods with spherical-like ends. 
Nonetheless the conclusion that large tori are energetically and geometrically favoured 
is a useful one, as the optimum micelle has a finite aggregation number. This is 
important. Otherwise, as we shall show, cylinders would continue to grow without 
bound, a prediction often in conflict with observation. 

For completeness, we now give the equations which describe the transition micellar 
shapes from spheres to cylinders, and examine some of their properties. 



(d) THE SPHERE-CYLINDER TRANSITION SHAPES 

The simplest transition shape consistent with all three packing criteria are ellipses 
of revolution, depicted in fig- 7. The ellipse is centred at a distance r from the origin 
and has major and minor axes equal to c and b. At the point P, where the curvature 
will be a maximum, the principal radii of curvature are 

Ri = r+b 9 R 2 = c 2 /b. • (5.6) 

For this shape the following conditions must apply : 

N = V/v = AJa 0 [criterion (2)] (5.7) 

/ = R 2 -» c*fb < 4, (5.8) 

where / the length of the hydrocarbon tail of an amphiphile at P. The assumption 
/ = R 2 is not obvious. Intuitively, the reason for this choice can be seen in fig. 5 and 
6. A truncated pyramid (or cone) would lead to a micelle with a region behind the 
hydrocarbon tail end of the amphiphile which cannot be filled up, Le., there results a 
micelle with a hole in it. Setting R 2 - / gives instead a wedge-shaped region (fig. 6) 
which avoids this difficulty. The packing equation [eqn (5.4)] determines another 
relation of the form 

R x =(r+fe)= 1 



3[l-2t>/fl 0 Q 
or 

b " 3[l-2p/a o rf { } 

The transition shapes may be divided into a globular regime and a toroidal regime. 

(e) GLOBULAR SHAPES (f < b) 

For this shape we find the following expressions for the total volume Kand surface 
area A of the hydrocarbon core : 

V - 27rc[r 2 sin 9 0 +rb(9 0 +i sin 20 o )+6 2 (sin 0 O -J sin 3 0 O )] (5.10) 

A = 4rcc|^^ (5.11) 

where 

cos 0 O = -(rjb). (5.12) 

The expression for the volume Vis exact, while that for the area A is approximate to 
first order in [l-(&/c) 2 ], which is entirely adequate for our problem since (bjc) is 
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always close to 1. Applying the condition given by eqn (5.6), (5.10) and (5.11), we 
obtain the following equation: 



Ww 



[i sin 0 o (2+cos 2 0 Q )-e 0 . cos 0,,] 



(5.13) 



Thus for a given value of i>/a 0 /, eqn (5.9), (5.12) and (5.13) may be solved to find the 
unique value of (b/c) from which all the other parameters such as b, c, r and N may be 
readily calculated. It can be seen that N/(4iil 3 /3v) is a function of v/a 0 l. 



(J) TOROIDAL SHAPES (r > b) 

For this shape the total volume V and area A of the hydrocarbon core is given by : 

V - libber (5.14) 

A = 4ji 2 ct[1 -i[l - (6/c) s D, (5.15) 

where again the expression for the area A is to first order in [1 -(6/c) 2 ]. For this 
shape (6/c) is readily obtained and is given by 

3vfa 0 l 



(blef 



(5.16) 



from which the values of b, c, r and N may be readily calculated. 

For both the globular and toroidal shapes the values of N have been calculated 
numerically. For fixed v and a 0 , N is a decreasing function of /, so that its minimum 
value M occurs when / = /„ i.e., the maximum value of /. Since we require the 
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Fig. 8.— Plot of optimal aggregation number M against (vla 0 k) for transition shapes of the ellipsoid 
f revolution model. Me = AvI */3i> is the largest possible aggregation number for spheres consistent 
with packing. Lower curve corresponds to spheres with no packing constraints. 
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shape with the minimum aggregation number we take / = / c hereafter. The variation 
of Mj(4rdll3v) with vja 0 l e is plotted in fig. 8. The range of values vfaok from J to ^ 
correspond to the transition from the purely spherical to the purely cylindrical shape. 
The results show that as v/a 0 l e increases above i the spherical shape becomes pro- 
gressively more globular (see fig. 7), the two radii of the ellipse become progressively 
smaller than / c while the ratio bfc decreases below L The transition from the 
globular to the toroidal shape occurs when v/a 0 l c has reached 0.383 (or about 15 % 
higher than its value of i when the spherical shape is no longer viable). At this 
point b/c has reached its minimum value of 0.74, and the micelle contains about twice 
the number of amphiphiles in a spherical micelles of the same value of v and / c . 
As vfa 0 l e further decreases the toroid grows outward ; both b and c increase towards 
/ c and the ellipse approaches a circle once more. At t?/a 0 / c = J, the critical condition 
for cylinders is reached ; the ellipse is now a circle and both b and c equal / c . For 
each of these shapes it may be verified that no point within the micelle is farther from 
the surface than / c , and that the total curvature at the outermost regions is greater 
than at any other part, so that the packing equation is everywhere satisfied. 

As mentioned earlier, the above model does not depend on the way v and / c are 
calculated for an amphiphile. This is a separate matter, and has been analysed by 
Tanford, 1 * 29 who used the X-ray data of Reiss-Husson and Luzzati 38 to obtain 
the following relations 

v « (27.4 + 26,9 n) A 3 per hydrocarbon chain) ^ 
/ c = (1.5+1.265 n) A per hydrocarbon chain j 

where n is some value close to but smaller than the number of carbon atoms per 
chain. Tanford 29 used the above relations to interpret some experimental results of 
Swarbrick and Daruwala 39 who measured the aggregation numbers of n-alkyl 
betaines in water as a function of aikyl chain length from C 8 to C 15 . He calculated 
the surface area per amphiphile on the basis of an ellipsoidal model and found, 
somewhat unexpectedly, that this area fell considerably with increasing chain length. 

Table 1.— Area per head group calculated for n-ALKYL betaines using present model, 
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Experimental values of obtained from ref. (38). 



We have repeated these calculations based on our ellipse of revolution model (see 
table 1), and find that the surface area is now practically invariant, Le. 9 it is almost 
insensitive to chain length, as expected. We may note, too, that in Tanford's calcu- 
lations, the surface areas were calculated at a distance of about 3 A from the surface 
of the hydrophobic core, whereas in our calculations the surface areas were calculated 
at the surface of the hydrophobic core. This is consistent with our earlier conclusion 
(section 3) that it is at this interface that the surface area per amphiphile a Q should be 



J. N; ISRAELACHV1LI, D. J. MITCHELL AND B. W. NINHAM 1545 



determined, and that it is this area that should remain fairly constant with changes in 
curvatures and/or alkyl chain length. Our results for the n-alkyl betaines bear out 
this conclusion and also lend support to our micelle shape. 

6. NON-SPHERICAL MICELLES: THERMODYNAMIC 
CONSEQUENCES OF PACKING 

The thermodynamics of spherical micelle formation as discussed in section 4 
provides an encouraging preliminary test of our model. We now wish to carry out 
a more refined comparison between experiment and theory, modified to include 
packing constraints. The analysis of this section falls into two parts. (1) A discussion 
of globular micelles. Here the model gives even better quantitative agreement with 
observation. (2) A discussion of the behaviour of very large aggregates. Again we 
compare with experiment 

(a) GLOBULAR MICELLES 

As discussed in section 5, once aggregation numbers exceed a certain limit, it is 
no longer possible to pack amphiphiles into spheres. If the aggregation number is 
not too large, we can expect micelles to take a globular, near spherical, shape and 
now consider corrections to the theory which arise from this non-sphericity. 

(i) CHANGES IN SALT CONCENTRATION 

We proceed as in section 4, where it was shown that because the mean aggregation 
number N is close to the (energetically) 44 optimal *' aggregation number M t the area 
per lipid is also close to a 0 ; the consequent error in assuming a = a 0 is not large. 
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For globular micelles it would be extremely difficult to attempt a complete 
thermodynamic analysis. Therefore, we shall assume that the micelles are narrowly 
dispersed and take on the optimal size and shape discussed in section 5, Le. 9 the 
shape and size of the smallest micelle which has a head group area a 0 . The only 
consequence then of packing constraints is that the relation between a 0 and M is 
modified. This relation is exhibited in fig. 8. Clearly, the value one obtains for a 0 
depends on what one assumes to be the largest aggregation number possible for a 
spherical micelle. For SDS we have taken this cut-off size (AT c ) to be either one of 
the two smaller values 57 or 64. Were we to take M c > 123 the graph of fig. 3 [in 
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which Oo/(^)*(3»)* is taken to be precisely M~*] would be reproduced. The results 
for M c ~ 57, 64 are tabulated in table 2 together with values of K u K 2 and y which 
best fit the experimental data exhibited in fig. 9. K u K 2 are defined by 



The best fit for M c = 64 corresponds to an interfacial energy at the micelle-water 
interface of y « 52 erg cm"" 2 in agreement with the measured value of y t at the bulk 
oil-water interface. The interfacial energies of liquid hydrocarbon-water interfaces 1 3 
vary from about 50 to 54 erg cnr 2 (from section 4 with the spherical approximation, 
the best fit y was 37 erg cnr a ). With a' ~ 20 A, the value of K 2 corresponds to a 
hydrophobic energy of 13 300 cal moH. Curvature corrections (see below) reduce 
this to about 12 000 cal moH close to the expected value deduced in section 4. 



Fig. 9.— Variation of ln(cme/55.5) with ao/Ofap*)* for SDS in Nad Data listed in table 2. 
O, Afc = 57 ; #, M c = 64 ; - - M c > 123 (this line ignores packing and is equivalent to fig. 3). 

(H) EFFECT OF CHANGES IN CHAIN LENGTH 

Our model can be subjected to a further test From eqn (6.1) for In (cmc) we 
find that the only term which depends on n is the hydrophobic term g—g'. It then 
follows that In (cmc) must be a linear function of n. This is found to be so experi- 
mentally. The experimental results for In (cmc) for non-ionic surfactants or ionic 
surfactants at 0.5 mol dm ~ 3 salt, correspond to a hydrophobic free energy 4 of about 
720 cal mol- 1 (CH 2 group)" 1 . The expected hydrophobic energy as deduced from 
solubility data for hydrocarbons in both water and in hydrocarbon is about 820 cal 
mob 1 (CH 2 group)- 1 . 4 This discrepancy can be resolved if we specialise our model 
to include the curvature effects discussed in section 3. Unfortunately there is some 
uncertainty as to the correct value of D which makes it difficult to give a proper 
quantitative account of curvature corrections, but we can attempt a rough accounting 
as follows. First, let us determine a range of values for D. We have from eqn 
(4.3) a 0 = y/ln^D/ey and from table 2 we find 



In (cmc) = Pyfoo-aO+tf-^W 



(6.1) 




.19 .20 .21 .22 .23 .24 .25 .26 .27 



0.22 < 



(4*)*(3p)* 



< 0.26, 



(6.2) 
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where for SDS *• 11 v « 350 A 3 . Whence D is in the range * 

l(T 9 <-<1.4xlO- 9 cm (6.3) 

over the range of salt concentrations 0 to 0.3 mol dm" 3 . Taking e = 40 we have 
4 < D < 5.6 A, and for e = 80, 8 < D < 11 A. The first choice seems more 
reasonable as the choice c = 40 corresponds more nearly with expectations on the 
dielectric constant near the oil-water interface. Now, assuming an approximately 
near spherical shape, the curvature correction to be added to eqn (6.1) is 

lal-^— -l)lkT= - , (6.4) 

\l+D]R Jl kT(R+D) % v ' 

where R ~ 3v/a 0 is the mean radius of the globule. This expression assumes that 
a a 0 . The parameters?, a 0 , D are independent of n> and R is a linear function of n. 
As a typical example, take a Q = 60 A 2 corresponding to M G = 64, y — 50 erg cm~ 2 , 
and D ~ 5 A. Then we obtain for various values of n using v = 27.4 + 26.9 n A 3 
[eqn (5.19)], the results in table 3. It can be seen from the table that curvature 
corrections are in the right direction in that they reduce the hydrophobic energy and 
are qualitatively correct whatever choice is made for Z>. 

Table 3 
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1.82 
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0.14 


66 


84 
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16.1 


1.72 


2.78 


















0.11 


0.15 


66 


90 


13 


17.5 


1.61 


2.63 


















0.89 


0.12 


54 


72 



For alkyi sulphates in the absence of salt both D and a 0 will be larger, leading to 
larger curvature effects. The observed discrepancy here however is about 400 cal 
moh* (CU 2 group)- 1 , 4 and this is too large a discrepancy to be accounted for by the 
simple capacitance model. At very low salt concentrations with ionic surfactants we 
expect highly non-linear behaviour for the repulsive electrostatic terms, so that this 
simple model can be expected to break down. 

(HQ TEMPERATURE EFFECTS 

We turn now to the temperature dependence of the cmc. Eqn (6.1) gives little 
indication of how the cmc should vary with temperature. The observed cmcs all 
have a minimum between 20-30°C, for sodium n-dodecyl, n-decyl and 2-decyl 
sulphates, and for n-dodecyl trimethylammonium bromide. 40 " 42 While the tempera- 
ture dependence of the term 2ya 0 /kT could be estimated from data on the surface 
tension of water and decreases with temperature, we have no information on the 
temperature dependence of the hydrophobic term g— g\ However, evidence 4 on 
the solubilities and enthalpy of solution of the sequence ethane, propane, butane, 
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together with corresponding data for the saturated alcohols from ethanol to n- 
pentanol, indicates that d$T[{g-g')lkT] increases as the number n of CH 2 groups 
increases. Since lya^ is independent of n, eqn (6.1) predicts that the slope of the 
cmc against T curve should increase as n increases. This qualitative prediction is in 
accord with Flockhart's results 40 on sodium n-alkyl sulphate micelles for n = 10, 
12, 14. 

One further prediction from eqn (6.1) is that the slope of the cmc against T curve 
should increase as the ionic strength increases. This follows because an increase in ionic 
strength causes a decrease in a Q while other parameters in eqn (6.1) remain constant 
Such an effect has been observed by Matijevii and Pethica 43 with SDS micelles in 
0.01-0.20 mol dm~ 3 NaCl solutions. 



When a 0 is small so that the optimal size is much larger than the largest possible 
sphere, we have argued in section 5 that the amphiphiles would pack into large 
toroids if the optimal si2e were indeed attained However, as can be seen from 
fig« 8, aggregates much smaller than the optimal may have an area per molecule, a, 
only marginally larger than a 0 and, consequently, self-energies $ only marginally 
higher. It follows from the arguments of section 2 that in this case the system will 
be polydisperse and we may expect to find aggregates somewhat smaller than the 
optimal size. This is in contrast to the situation with globular micelles and with 
lipid vesicles to be discussed in the next section. The problem of determining the 
shape of the micelles is much more difficult For micelles of a given aggregation 
number N 9 and, therefore, of a given volume Nv, and a given total area A, the 
minimum $ can be shown to correspond to uniform area a - A/N, but only provided 
this shape is allowed by packing. The minimum $ for the given aggregation number 
N would be determined by the smallest value of a. The functional relationship 
between this minimum area a and the given aggregation number N is the same as that 
between a 0 and the minimum aggregation number N illustrated in fig. 8. Use of this 
form leads to results which do not agree with the observed behaviour of large 
aggregates. However, for some values of N packing restrictions may prohibit the 
formation of a micelle with uniform area a, and the micelle with minimum jtft may 
have a non-uniform area per molecule. A possible candidate is a cylindrical micelle 
with a = oq having globular ends with a > a 0 . Such a shape would have an average 
area per molecule very close to a 0 . Experimental evidence seems to indicate that 
micelles of large aggregation number formed from one-tail-amphiphiles are rod- 
shaped. 36, 37 Such a micelle would have a of the form 



where f& is the self-energy per molecule of an infinite cylinder = 2ya 0 . The term 
oJcTJN arises because the molecules in the globular ends have a larger area and the 
number of these is independent of N. As was shown in section 2, this form for $ 
leads to very polydispersed systems whose mean aggregation number increases rapidly 
with amphiphile concentration, provided a is large. It is experimentally observed 
that the mean aggregation number of these (apparently rod-shaped) micelles does 
change rapidly with concentration, eventually levelling out to a fixed number. 4 

A model which might be expected to give a rough estimate of a is shown in fig. 10. 
Consider a cylinder of radius r = 2vja 0 with spherical ends of radii 4- Then the 
molecules in the hemispheres on the ends have area a = 3i?// c and all other molecules 



(6) ROD-SHAPED MICELLES 




(6.5) 
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in our model have a = a 0 . Then if tj « 4*/? /a denotes the number of molecules in a 
spherical micelle of radius / c we find that $ is given by 



with 



This gives 



with 



(6.7) 



(6.8) 



x « flo /a = * 0 ' c /3i; (6.9) 

so that x lies between 1 (corresponding to the smallest a 0 attainable by spheres) and 
$ (corresponding to the smallest a Q attainable by infinite cylinders). 

In table 4 we tabulate values of a for various values of x, assuming different values 
for / c with y = 50 erg cm- 2 . 
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Fig, 10.— Cylindrical micelle with globular ends. Aggregation number can be less than that of the 

optimal torus. 
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TABLE 4 

— 15 A , Jo-25A 

* a a 

1 0 0 

0.95 0.9 2.4 

0.90 3.5 10.0 

0.85 7.9 21.8 

0.80 14.0 38.8 

0.75 21.8 61 

0.70 31.4 — 

0.60 38.8 — 

Table5 



* -. 



* J> 

3 



1550 



ASSEMBLIES OF HYDROCARBON AMPHIPHILES 



Now it has already been shown in section 2 [of. eqn (2.16)] that for S > cmc and 
a > 1, N = 2^/Se* (as was obtained by Mukerjee, 44 by a different approach). 
In table 5 we tabulate N for various values of a at different concentrations in mole 
fraction units. 

Some interesting features of this model deserve comment The values of a which 
best fit available experimental data were found by Mukerjee 44 to lie in the range 
25-35. For the amphiphiles concerned / c ~ 20-25 A. It can be seen from tables 4 
and 5 that the model used here does in fact produce values of a which straddle this 
range. An increase in salt concentration decreases a 0 and therefore x, so that N 
increases rapidly with increase in salt concentration. This is in keeping with 
experimental observation. 36 Also an increase in chain length increases / c and 
\ consequently increases a. Therefore amphiphiles with larger chain lengths will tend 
^ to have micelles with larger aggregation numbers which increase faster with 
concentration 4 

There is a further interesting consequence. As the lipid concentration is increased 
and the aggregation number is, as a consequence, increased, the " optimal size " 
will be approached. After this value $ will remain constant and the micelles will 
cease to grow rapidly. This seems to be in keeping with experimental observation 
although there is some disagreement about interpretation of the results. 4 * 44 

(c) BILAYERS 

It is of interest to note here that the behaviour of lipid aggregates for which the 
optimal shape is a spherical vesicle or an infinite bilayer will be quite different from 
that of rod-shaped micelles. One might anticipate by analogy with rods that finite 
disc-like bilayers could form with cylindrical rims. For such a system the free energy 
would take the form /$ = p% + afcT/NK From section 2, for large a, the system will 
assemble spontaneously into infinite bilayers. Repetition of an analysis similar to 
that for rods leads to values of a sufficiently large so that this process always occurs. 

7. THE STRUCTURE AND STABILITY OF 
ONE-COMPONENT BILAYERS 

The preceding section has established that the predictions of our model accord 
both qualitatively and quantitatively with available observations on globular and 
rod-like micelles. We can, therefore, now turn with some degree of confidence to an 
area as yet imperfectly explored, but whose importance is not in dispute. We shall 
apply packing criteria to a consideration of those amphiphiles which coalesce naturally 
into vesicles or planar bilayers. Attention will be confined to biological phospho- 
lipids : these have anionic or zwitterionic head groups; they commonly contain two 
hydrocarbon chains with 16-18 carbons per chain, and display some unsaturation. 
It is this unsaturation that ensures that the hydrocarbon chains are above their 
" melting temperature i.e., that they are in a fluid-like state, above 0°C 4 ' 7 * 45 
Because these molecules have two hydrocarbon chains, they tend to have double the 
volume of the single chain molecules previously considered, for which vfa c l e £ \ ; 
thus for phospholipid molecules v[a 0 ! c > and globular and cylindrical micelles are 
prohibited. 

We first analyse the free energy of a one-component spherical vesicle bilayer 
(fig. 11), and investigate its stability in the absence of packing restrictions. It should 
be emphasised that without packing no sensible results emerge unless we include 
curvature corrections. Subsequently, when packing is built into the model, curvature 
corrections become of secondary importance, so that the principal results do not 
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depend on a detailed model for the repulsive stabilising forces. In accordance with 
the terminology of the preceding sections, the free energy per amphiphile in a spherical 
vesicle of outer and inner radii JJ, and R 2 is given by 

where and n 2 are the number of arophiphiles in the outer and inner layers of the 
bilayer, e is the charge per polar head group, a is the dielectric constant of the head 
group region, and D is the effective capacitor thickness of the head group. The 
interfacial energy of the hydrocarbon-water interface is given by 7, and the outer and 
inner radii, R t and R 2 , are referred to these two hydrocarbon-water interfaces. 
Eqn (7.1a) is written in the more convenient form eqn (7.16) in which the elusive 
parameters e f D and a are replaced by y and the optimum surface area per amphiphile 
in the bilayer a 0 , as given by eqn (4.3). 




Ro. 11.— Spherical bilayer vesicle of hydrocarbon thickness /. 

The aggregation number N and hydrocarbon thickness t of the vesicles are given by 
N = n 1 +n 2 — 4n(Rl -Rl)/v (7.2) 

and 

t = R 1 -R 2 , (7.3) 

where v is the hydrocarbon volume per amphiphile, it being assumed that, to a first 
approximation, v is constant and independent of vesicle size. If a x and a 2 are the 
surface areas per amphiphile on the outer and inner hydrocarbon-water interfaces, 
we may write 

a, - 4nRUn t (7.4) 

and 

a 2 = 4zRl/n 2 . (7.5) 

Eqn (7.1)-(7.5) allow us to determine the minimum free energy configuration for a 
particular value of NorRx. This is done in Appendix B where, in the absence of any 
packing restrictions, we obtain the following results : 
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j&min) « 2a A 



a i ® <*o\ 1 



a 2 « a, 




(7^6) 



(7.9) 



(7.7) 



(7.8) 



where a 0 is the surface area per amphiphile in the planar bilayer (i,e., when N = oo or 
iJj = oo). From the above results we may conclude the following : (1) the minimum 
free energy per amphiphile is slightly lower than that in a bilayer, for which /iftOain) = 
2a 0 y. Since a spherical vesicle must have a much lower aggregation number N than 
a planar bilayer we conclude that spherical vesicles are thermodynamically favoured 
over planar bilayers and that, in the absence of any packing constraints, the vesicles 
will have very small radii. Indeed, if there were no packing constraints the vesicles 
would shrink to such a small size that they would actually form into micelles. As 
will be shown, packing constraints prevent vesicles from shrinking much below a 
certain critical packing radius. (2) The optimum area per amphiphile on the outer 
surface a t is found to be smaller than the optimum area for a bilayer a 0 , while the 
inner area a 2 is found to be larger. (3) The hydrocarbon thickness t is slightly greater 
than that of a bilayer, given by t = 2v/a 0 . For example, if t = 30 A, R 1 = 100 A, 
the value of t is increased by about L5 %, 

While some of the above results are modified when packing constraints are 
introduced, the important conclusion that remains unchanged is that spherical 
vesicles are thermodynamically favoured over planar bilayers. The effect of packing 
is only to place limits on their size. We may note here that while sonication is generally 
employed in the laboratory for their formation this is not essential : Batzri and Korn 46 
have reported the formation of vesicles without sonication. In addition, vesicles, 
once formed, appear to be stable for many days. 46 ' 50 We return to this matter later. 

(a) THE EFFECT OF PACKING CONSTRAINTS ON BILAYER STRUCTURES 

In section 5, a geometric "packing equation" was derived that relates the 
amphiphile surface area a, hydrocarbon volume v, and hydrocarbon thickness / to 
the curvature at the amphiphile surface. For a spherical vesicle of outer radius Ri 
the packing equation, eqn (5.4), for the outer layer amphiphiles is 



It will be shown subsequently that the thermodynamically favoured vesicles have their 
hydrocarbon chains maximally extended. Then / = / c , and R t takes the minimum 
possible value for a given a t . This value will be called the critical packing radius R\. 
To a first approximation the outer amphiphile surface area a x may be put equal to a 0 » 
and the value of R] may then be readily obtained from eqn (7.1 1). As an example of 




(7.10) 



so that 



6(1-1^0 



(7.11) 
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eqn (7.11) we shall use some accurately determined data 51 » 52 for egg phosphatidyl- 
choline (egg PC) for which v = 1063 A 3 , a 0 = 71.7 A 2 , and we obtain the following 
critical vesicle radii as a function of / c . R\{! c = 19 A) « 80 A, R\(l G = 17.5 A) « 
108 A, = 16A)« 212A. 

Before we proceed with a thermodynamic analysis it is worth considering the 
qualitative consequence of geometric packing in order to understand why double 
chained arophiphiles aggregate so differently from single chained amphiphiles. In 
section 5 we found, as a rule of thumb, that those amphiphiles for which p/a 0 4 lies 
between } and | form into micelles whose shapes vary gradually from spherical to 
cylindrical. For values of v/a 0 ! e between i and 1 we might expect a gradual variation 
from a cylindrical micelle to a planar bilayer. However, there is no way that a 
cylinder may be gradually distorted into a planar bilayer consistent with the geometric 
packing criteria. On the other hand, once the critical condition for cylinders is 
reached, i.e., once i>M)/ c = i» it may he readily verified that a very small spherical 
vesicle may already be packed with an outer radius Rx equal to (1 + 1 / ^/ 3)/ c « 
1.6 l c . Thus we do not need to look for a gradual transition from cylindrical micelles 
to bilayers. At some value of v}a 0 l c close to \ the spherical vesicle will become 
viable as regards packing, and since its aggregation number is bound to be much 
smaller than that of a long cylinder it will take over as the most favoured structure. 
As v/a 0 l c increases above \ the vesicle will grow, becoming a planar bilayer when 
v/a 0 l e reaches 1. 

Tanford *• 29 has concluded that almost all of the commonly studied amphiphiles 
with one hydrocarbon chain cannot pack into spherical micelles. Such amphiphiles 
invariably form into small globular micelles or long rod-like micelles. Thus their 
values of v/a 0 l c lie in the range J to If the only effect of doubling the number of 
hydrocarbon chains is to double v[a 0 l cf we may conclude that the values of v[a 0 l c for 
diacyl chained phospholipids should lie in the range f to 1.* By this criterion, using 
eqn (7.1 1), such phospholipids should form either planar bilayers or spherical vesicles 
of radii i? x > ~ 3/ c . 

(b) THERMODYNAMIC CONSEQUENCES OF PACKING CONSTRAINTS 

For vesicles that are larger than the critical radius R u there are no packing 
restrictions on any of the amphiphiles and the minimum free energy configuration is 
given by eqn (7.5>(7.9). Such vesicles, however, are not thermodynamically 
favoured over smaller vesicles which have lower aggregation numbers N. The 
vesicles will therefore tend to shrink ; but once their outer radii become lower than 
the critical packing radius R\, the surface areas of the outer amphiphiles a x must 
now increase above the optimal area [given by eqn (7.7)) and this increase 
in a x leads to an increased free energy per amphiphile as N or R y falls. We 
must, therefore, analyse how the free energy per amphiphile varies with N when 
packing constraints are included* To obtain exact solutions to the initial eqn. 
(7.1)^(7.4) is difficult when packing constraints, given by eqn (7.10), are introduced* 
But since we are concerned with vesicles where the effective capacitor length D is 
much smaller than the outer and inner radii R x and R 2 approximate solutions may 
be obtained by ignoring curvature effects, which are now small. This is done in 
Appendix B where it is first shown that the optimum surface area per amphiphile on 

* As already discussed, the assumption that the surface area a 0 is not much affected by doubling 
the number of chains is supported by the observation that these areas are independent of chain length, 
and have much the same values for both single chained and double chained amphiphiles. This is 
also expected on theoretical grounds for hydrocarbon chains in the fluid state. 
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the inner surface a 2 is unaffected by the packing constraints, which, as is to be 
intuitively expected, only affects the outer layer amphiphiles. The inner layer | 
amphiphiles are always in a favourable packing state; their hydrocarbon chains } 
simply fill up the inner layer volume with their inner surface areas a 2 remaining at 
their optimal values, close to a 0 . In Appendix B it is further shown that ignoring 
curvature effects the mean free energy per amphiphile is given by 

4ji/ 2 v 

tivtoon-fd-RtlRl) 2 (7.12) 

when R t is smaller than the critical packing radius R\ given by eqn (7.1 1). For 
R x > R* we have 

$*2a 0 y. (7.13) 

To obtain the vesicle size distribution we substitute the above expressions into 
eqn (2.8). For R x < R\ we obtain the following distribution law for the concentra- 
tion of vesicles X Rl of radius R x : 

(x » \w* 
-jZ) exp {-4*fr(l -RjRtflkT], (7.14) 

where NfN c is the ratio of the aggregation numbers of vesicles with radii R x and HJ, 
and is given approximately by 

NIN< - {^+(J? 1 -0 2 }/{(*!) 2 +(R!-0 2 }. (7-15) 
When Ri > R? x the distribution function reduces to 

Eqn (7.14)-(7.16) give the concentration of vesicles as a function of R x . As an 
application of this distribution we shall use the values previously obtained for egg 
phosphatidylcholine. Fig. 12 shows how the egg PC vesicle concentration varies 
with outer radius Rx for two values of l e and their corresponding values of R\ 
determined by eqn (7.11), i.e,, / c = 17.5 A, R\ = 108.3 A, and / c - 16.0 A, R\ « 
212.5 A. We have assumed that y - 50 erg cnr 2 , / = 29.6 A," T « 20°C and 
have plotted the distributions for (X^JN 0 ) « 10~ 8 and (X^JN') = 1(H 4 . The 
value of (XjfJN c ) = 10~ 8 corresponds roughly to a concentration of 1 mg phospho- 
lipid per cm 3 water. 

A number of general conclusions for egg PC may be drawn from the distribution. 

(1) The concentration of vesicles has a near-Gaussian profile which peaks at an 
outer radius R t close to, but slightly smaller than, the critical packing radius R\> 
At the peak, the surface area per amphiphile is slightly greater (by ~1 %) than that 
of the bilayer a 0 , though the difference is diminished as the total amphiphile concentra- 
tion increases. 

(2) The profile is given approximately by exp {-(R x -R[) 2 f2c 2 } f where the 
standard deviation c is given by 

Thus the standard deviation a is found to be proportional to the vesicle size R\, and 
is approximately equal to 3 % of R\, Further, the theoretical value of a is fairly 
insensitive to the value chosen for the hydrocarbon-water interfacial energy y. 
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(3) The mean outer radius decreases as the lipid concentration is lowered, but the 
magnitude of the shift is very small over a wide range of concentrations. A million- 
fold change in concentration shifts the peak by only a few Angstrom units. 

The implications so far are that egg PC vesicles should be fairly homogeneous in 
size, and that this size is effectively independent of the total amphiphile concentration 
even down to extremely low concentrations. In general, the factors that are expected 
to affect vesicle size are changes in the optimal surface area a 0 and the ratio vf! G . 
For anionic amphiphiles the surface area a 0 may be expected to be changed by varia- 
tions in pH, ionic strength, [Ca 2+ ], e/c., 53 while zwitterions may be expected to be less 
sensitive to such changes. On the other hand, v/l c may be expected to vary with 
the degree and type (cis, trans) of unsaturation. 



4* 17.5 A *>16.0A 




vesicle radius, Rijk 

Fig. 12.— Concentration of phospholipids Xn in spherical vesicles of outer radii R x . The curves 
are based on eqn (7.14X7.1Q and are plotted for y = 50 erg cm -2 , t = 29.6 A (corresponding to 
egg PQ, for two values of / c . The vesicle concentration (Ajfi /tf c ) is in mole fraction units (55.5 mol 

dm- 3 ). 

We now return to fig. 12 in order to test its predictions for egg PC vesicles. These 
vesicles have been found experimentally 54 to have a total outer radius of 105+4 A. 
(A knowledge of the exact outer radius is not essential in the foregoing analysis.) 
Since the outer radius R t is measured at the hydrocarbon-water interface, a total 
radius of 105 A corresponds to a value of R t closer to 100 A or even less (depending 
on the length of the polar head group). Fig. 12 shows that for the vesicle size 
distribution to peak at R t « 100 A requires the maximally extended hydrocarbon 
length to be / c « 17.5 A. The theoretical results also predict that the vesicle 
diameters will have a standard deviation of ~1 A, which may be compared with the 
experimental value of ±8 A, 54 and that the aggregation number (at R t = 100 A) is 
2600 as also observed. 54 Further, the ratio of the number of egg PC molecules on 
the outer and inner vesicle layers should be given roughly by [R\}{R\— t)] 2 , where t 
is the hydrocarbon thickness and not that of the whole vesicle bilayer. Putting 
/ = 29.6 A yields an outer PC/inner PC ratio of (100/70.4) 2 = 2.0. If instead of 
29.6 A a value of 40 A were chosen for t (as is commonly done) we should have 
obtained a ratio of (105/65) 2 = 2.6, while for t = 45 A the ratio rises to 3.1. We 
may note that experimental determinations of this ratio vary greatly and depend on 
the method of measurement. The reported values to date are: 1.6, 48 1.7, 55 1.85- 
1.95, 5 *- 59 2.0, 60 2.2, 61 2.3, 49 - 50 2.6 9 (the average value is about 2.0). 

As regards the outer and inner surface areas of the amphiphiles, the results 
indicate that these should be about the same, and close (within ~ 1 %) to the area 
in lamellar bilayers, as observed. 9, 62 
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The value obtained for the maximally extended hydrocarbon length of 17.5 A is 
more difficult to compare with experimental values. In section 5 we found that 
80-100 % of the fully extended hydrocarbon chain appears to be a reasonable 
estimate of / c of saturated chains (see also Tanford 4 > 29 ). Since egg PC is mainly 
16:0, 18:0 and 18: 1, 63 we obtain a theoretical estimate for l e of (0.8-1.0) x 
[1.5+(1.265 x 16)] « 17.4-19.4 A. Alternatively, the fully extended length may be 
obtained from the limiting surface area at which monolayers collapse. For egg PC 
this occurs at a limiting area 64 » 65 of 55-62 A 2 corresponding to a fully extended 
hydrocarbon region (assuming v = 1063 A 3 ) of / c = 17-19 A. (We note that the 
value of / c at collapse may be somewhat greater than in the vesicle bilayer since the 
collapse surface area is appreciably less than the vesicle surface area. The difference 
is about 10 A 2 and implies that at collapse there will be some additional hydrocarbon 
chain available to the interior region.) 

The above conclusions appear to account fairly well for the experimentally 
observed behaviour of egg PC vesicles. The qualitative features of the results, 
however, should be generally applicable to all one-component vesicles composed of 
amphiphiles whose hydrocarbon chains are in a liquid-like state. 

A more detailed analysis of the problem, that includes the electrostatic effects of 
curvature, is difficult and probably unnecessary. In section 4 the effects of curvature 
on micelle size were analysed and it was found that, depending on the value of the 
dielectric constant e, the micelle distribution was shifted to slightly larger sizes (fig. 2). 
In the case of vesicles a similar result is expected. The electrostatic effect of curvature 
will lower the optimum area per outer amphiphile below a 0 as R t decreases [see 
eqn (7.7)]. This will shift the critical packing radius R\ to larger values and thereby 
shift the distribution curve to larger radii. On the other hand, the effect of curvature 
should lower the minimum energy per amphiphile as R x falls [see eqn (7.6)], and 
this will shift the distribution curve to smaller radii. These two effects will partly 
cancel out and the net shift could be in either direction. If, as expected, the net shift 
is to larger radii the implication is that / c is slightly greater than 17.5 A (from fig. 12 
we see that a small change in / c causes a large shift in the distribution curve). 

Finally, the above analysis allows us to calculate a value for the Young's modulus 
of a planar bilayer under compression or tension. Rewriting the free energy per 
lipid, in terms of the volume v and hydrocarbon thickness / (a = v/t) f the energy per 
two lipids on either side of a bilayer is = 2y(a+a$/a) ■» 2y(v( t+talfv). Thus 
the force F on compressing the bilayer from its initial equilibrium thickness f 0 to 
(t 0 - At) is 

at \ t 2 v) 




«4yo 0 ^. (7.18) 
'o 

The Young's modulus of the bilayer is therefore 

Af/t 0 t 0 

i.e., Yt 0 fa 4y « 200 dyn cm" 1 



(7.19) 
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For a bilayer of total thickness 50 A the effective Young's modulus is of order 4x 
]0 8 dyn cm~ 2 when an interfacial energy ofy w 50 erg cm' 2 is taken. Experimentally, 
it is only possible to measure Yt 0 . Though experimental determinations of Y vary 
greatly, Requena et a/. 66 have recently argued that the black lipid membranes have 
a value of Y close to 2 x 10 8 dyn cm -2 . The theoretical value of Yt Q « 200 dyn cm" 1 
may also be compared with measured values in the range 51-357 dyn cm- 1 for red 
cell membranes. 67 

We note that the elasticity discussed above is only for planar bilayers under 
compression or tension and does not extend to the bending of bilayers. On the 
contrary, our analysis has shown that the bending of a bilayer is favoured down to 
the critical packing radius (assuming that the lipids can freely rearrange by lateral 
movement and/or flip-flop), and that " bending elasticity " sets in only for radii 
smaller than this critical value. The elasticity of a fluid bilayer is therefore seen to be 
profoundly different from that of a classical elastic plate or shell. 



(c) FURTHER ASPECTS OF BILAYER STRUCTURES 

In general we find that, when amphiphiles aggregate, entropy favours the smallest 
structures, but that packing restricts the sizes and shapes of these structures. Single 
chained amphiphiles form small globular micelles or long rod-like micelles, whereas 
double chained amphiphiles form bilayers. But there are other important differences 
between single and double chained amphiphiles that are of some biological significance. 
The first is to do with the large dispersity of micellar size and aggregation number, 
and the great sensitivity of micellar size to changes in the total amphiphile concentra- 
tion. By contrast, double chained amphiphiles form very homogeneous structures 
that are unaffected by the total amphiphile concentration. Further, while single 
chained amphiphiles have relatively high critical micelle concentrations, of order 
10" 3 mol dm- 3 , double chained amphiphiles such as biological phospholipids, have 
very low critical micelle concentrations, of order 10" 10 moldm -3 . 4 It is these 
properties that make phospholipids such excellent building blocks for the formation 
of homogeneous and stable bilayer membranes, unaffected by drastic changes in the 
surrounding amphiphile concentration. However, there is always a continuous 
exchange of lipids with other membranes that is largely energy independent, 68 * 69 and 
which may be important in membrane-membrane interactions. 70 

In the other extreme of very high amphiphile concentrations a very different 
picture emerges as bilayers (or micelles) now begin to interact with each other. 
We stress that our treatment applies only to dilute systems, where there are no 
interactions between the aggregated structures. Thermodynamically, this was 
expressed by putting the activity coeflScient equal to unity in eqn (2.1). If such inter- 
actions are present the expressions for the free energy in both the aggregated and 
dispersed states will be modified by additional terms due to long range van der 
Waals and electrostatic interactions. Such interactions become important at high 
concentrations and lead to the formation of ordered mesomorphic phases. 7 * 51, 52, 71 " 73 
Our present treatment clearly does not apply to these structures, though some of the 
observation on mesomorphic phases have a bearing on bilayer vesicles. One of 
these concerns the observation that phospholipids do not always readily form into 
vesicles, but rather into multilamellar bilayers or myelin figures. Sonication is often, 
but not always, necessary to break up and transform these bilayers into vesicles. 
But once formed, these vesicles are homogeneous and stable, and unaffected by the 
length of time and intensity of sonication. 54 It seems that as soon as bulk phospholipid 
is placed in aqueous solution the lipids initially form into the " neat mesomorphic 
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phase " (/>., bilayers) ; 7, * u s2 > 71 » 72 since as the lipids diffuse out of the bulk into 
the water the local lipid concentration is initially very high. Once formed, the 
bilayers cannot reform into vesicles without first having to break up. For this to 
occur, a certain activation energy must be overcome, and it is this energy that is 
furnished by the sonication process- But even then the total lipid concentration 
must be sufficiently low for the vesicles to be dilute enough not to interact significantly 
with each other. 

Thus whereas we conclude that spherical vesicles, if allowed by packing, are 
thermodynamically favoured over extended bilayers, we add the reservation that they 
need not necessarily form spontaneously from bilayers. 

It is worth examining how a vesicle might form spontaneously from a bilayer; 
this is shown in fig. 13. Our analysis suggests that the initial stages of the formation 
of a spherical bulge is energetically favourable [cf. eqn (7.6)], but the final stage of 
separation involves adhesion and fiision at the neck as the vesicle comes away from 
the bilayer. This often requires a certain amount of energy and brings us to the 
mechanism of membrane fusion which is beyond the scope of this paper. 



5 bilayer t 




Fig. 1 3. — Stages of vesicle and tubule formation from a one-component bilayer. Stages A and B are 
energetically favourable if allowed by packing, but require free lateral diffusion of lipids in the bilayer 
and water flow across the bilayer. Stage C involves fusion. Stage D is energetically unfavourable 
for a bilayer in the fluid-state. A mechanism of vesicle and alveoli formation from a membrane 
similar to that shown here has been found to occur in a variety of cell types and in micropino- 

cytosis. 74 . 75 

So far we have only looked at spherical or planar bilayers. More generally, the 
packing equation [eqn (5.4)1 shows that for a given value of a 0 , v and / cl the critical 
packing curvature, KW + W), is approximately constant (where R\ and R\ are 
the two outer radii of a curved bilayer at critical packing). Thus for a tubular vesicle, 
putting R* 2 =s oo, its critical radius F\ would be about half that of a spherical vesicle. 
Tubular vesicles, however, have a larger aggregation number than spherical vesicles 
so that they should not form and they have not been reported. In addition, we show 
in Appendix B that the free energy per amphiphile in a tubular bilayer is higher 
than that in a spherical or planar bilayer. Thus we conclude that tubular protrusions 
should not grow spontaneously out of bilayers (fig. 13), and that bilayers in a fluid- 
like state should not transform into tubules. For tubular membranes to exist they 
must either have a resilient cytoplasmic skeleton that gives them their structure, or the 
membranes must be rigid. Though this conclusion may not necessarily extend to 
multicomponent lipid bilayers or biological membranes, we have not found any 
evidence for unsupported tubular membranes that are in the fluid-like state. Thus 
tubular cell surface projections such as filopodia, flagella, cilia and microvilli all 
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appear to be supported by a cytoskeleton of rigid microfilaments or microtubules. 75 ' 82 
Those tubular structures that do not have a supporting skeleton appear to have rigid 
(non-fluid) membraneous walls, for example, microtubules, 82 invertebrate photo- 
receptor membranes, 83, 84 tubular mitochondrial cristae. 75 ' 85 ■ 86 

We now turn to the factors that affect vesicle properties. For anionic phospho- 
lipids, the optimum surface area a 0 may be expected to be larger at higher pH or 
lower ionic strength. 53 This will result in the formation of smaller vesicles, as 
observed for phosphatidylserine at different ionic strengths. 87 

A very , different matter is to change the properties of a solvent that already 
contains stable vesicles in suspension, for now any changes in pH, ionic strength, 
[Ca 2+ J, temperature, other lipid levels, e/c, will all induce strains on the already 
existent vesicles. These strains will be asymmetrical about the outer and inner faces 
of the vesicle bilayer, 88 and the way the vesicles will respond to these changes will 
depend on such factors as the flip-flop rate, the permeability of the bilayer to water 
aod ions, and the strength and energy of rupture of the vesicle. For example, on 
addition of Ca 2+ anionic vesicles are bound ultimately to rupture. 89 ' 90 The 
ruptured bilayers may then reform into larger vesicles or extended bilayers, 90 and it is 
probably this two-step process that the " fusion " of vesicles under the influence of 
Ca 2+ occurs. 91 On the other hand the extraction of Ca 2+ (by EDTA, for example) 
may cause existing vesicles to shrivel up. We have not attempted to analyse all the 
possibilities, though in many respects our conclusions are in agreement with our 
earlier model for the packing of lipids in bilayers, 92 and with the " bilayer-couple " 
hypothesis recently proposed by Sheetz and Singer. 93 

8. CONCLUSION 

The main aim of this paper has been geared towards the building and testing of a 
simple model for self-assembly of ionic and zwitterionic amphiphiles. This model 
exists on two levels. At a non-specific level, as underlined in earlier work by 
Tanford, 4 the principle of opposing forces, notions concerning geometric packing 
properly delineated, and thermodynamics, form a framework for a model of such 
fundamental and general application that its predictions cannot be seriously in 
dispute. It is reassuring that this indeed turns out to be so, both qualitatively and 
quantitatively. At the level of the simple model, in dealing with comparisons of 
theory and experiment on micellisation, the major new point of departure in our work 
has been the recognition of and an attempt to quantify the role of packing. When 
packing constraints are recognised a number of apparently disconnected facts begin 
to fall into place. Thus, for globular micelles, for the effects of salt changes on cmc, 
the best fit occurs for an interfacial tension at the oil-water interface close to 50 dyn 
cm- 1 , the known value. The rapid growth of rod-shaped micelles from globular 
micelles with changes of amphiphile concentration again emerges in a very natural 
and inevitable way once packing constraints are allowed. Further, not only do 
predicted parameters for the growth process agree with observation, but the packing 
model also predicts the observed levelling of mean aggregation number with increasing 
number per micelle, and the dramatic change in N with increasing chain length. 

At another and more restricted level, we have imposed much more stringent 
requirements, as curvature effects stem from a highly specific form for the repulsive 
electrostatic forces. It is again reassuring, although not a necessary requirement, 
that such a decoration removes at least in part two sources of difficulty. These are 
the annoying discrepancy of 100 cal moH (CR 2 group)- 1 between calculated and 
observed hydrophobic energies of alkane chains, and the apparent non-existence of 
tubular bilayer vesicles. 
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Especially pleasing is the circumstance that our model gives a unified framework 
which includes bilayers with no additional assumptions. Various physical properties 
of bilayer vesicles are readily accounted for. Even the Young's modulus of bilayers 
appears to be given correctly from elementary considerations. 

As far as simple modelling of self-assembly is concerned, the treatment of single 
component lipid molecules given here has probably been pushed as far as it can. The 
refinement of our theory of self-assembly requires a proper examination of Stern 
layers, consequences of deviations from liquid-like properties of hydrocarbon 
chains, 94 head group steric effects, specific ion adsorption and other effects. While 
such a more rigorous analysis would undoubtedly provide specific insights into the 
properties of particular molecules, it is doubtful if a more refined theory will provide 
a better overview. 

The simple model does allow an entry point into the study of self-assembly of 
multicomponent lipid systems, lateral phase separation (clustering), membrane 
asymmetry, and in particular how these relate to curvature through packing. These 
form a central class of problems in membrane biology. 

This work arose from some earlier unpublished ideas of Dr. L. R. White and 
one of us (J. N. I.), which were presented at the 49th National Colloid Symposium, 
Potsdam, N.Y., in 1975. 

APPENDIX A 
ELECTROSTATIC ENERGY OF ZWITTERIONIC ARRAYS 

We give here an example which illustrates how the capacitance form for the 
electrostatic interaction energy comes about for the special case of zwitterionic head 
groups. 

Consider a planar lattice of parallel M dipoles " oriented normal to the plane. 
The positive and negative charges of each " dipole " representing a zwitterionic head 
group are separated by a distance D, and the array is immersed in a medium of 
dielectric constant 6. Let a denote the area of an elementary cell of the lattice, Le. % 
the area per dipole. The electrostatic interaction energy per dipole is then 

where r denotes the position vector of a dipole relative to a fixed dipole and the sum 
ranges over the whole latti ce excludi ng the reference dipole (m = n = 0). For a 
square lattice rj= |r| = y/(m 2 +n 2 )a. We see immediately that u tl = e 2 /eDxa 
function of Df^Ja. In the limit v 9i reduces to the energy of a lattice of point 

dipoles of dipole moment eD. We are concerned with the opposite limiting case 
Djyja -+ oo, and first obtain the leading asymptotic term, which follows quite easily. 
Taking Mellin transforms 95 with respect to D, and inverting, we can estimate the 
integral representation 

~ 4^iJA3j r(i) "iW" • (A } 

where the contour of integration is a line parallel to the imaginary axis which satisfies 
-2 < Re j < 0. This integral representation can be substituted into eqn (Al) and 
the summation and integration interchanged provided we now choose y such that 
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~2< Rts< -1, since the sum Ir*- 1 converges when Rej< -1. Thus we 
obtain 

The asymptotic form of u tX as Z> -> oo is now given by the residue at the pole s *= - 1. 
To evaluate the residue we observe that the sum I r*~ l has the same singularity as the 
integral 

2n(a 



ijV. 



d 2 r 



(s+l) 

Evaluating the residue of the integrand of eqn (A3) at the pole s = — 1, we have 

2ne 2 D 

««i ~ • (A4) 

(D-oo) *<* 

The above argument does not depend on the assumption that the dipoles lie on a 
lattice since only the residue of the sum S r*- 1 at its pole at s = — 1 comes into the 
final answer, and this will be the same even if the dipoles are irregularly distributed so 
long as they have an average density \{a. Indeed, if we consider the correlation 
energy of a statistical distribution of dipoles we have 

where g(r) is the radial distribution function (which tends to unity as r oo). This 
expression leads to the same asymptotic form eqn (A4), and since the asymptotic 
result is linear in the coupling constant so, therefore, will be the free energy. 

To establish the range of validity of eqn (A4) we need further terms in a complete 
asymptotic expansion. This is a more difficult problem and to derive such expansions 
we return to eqn (A3). The sum can be carried out exactly 96 in terms of known 
transcendental functions to give 

where £(z) is the Riemann zeta function, and J?(z) is a related entire function 96 
defined for Re z > I as 

00 (—1)* 

The properties of fi(z) are enumerated in the important paper of Glasser. 96 Thus 

U * = -ip + 2^ra\X^) ^-^riJJ (A8) 

with -2 < Re^< -1. Translating the contour to the right we have poles at 
5 « - 1 with residue (—2) from the zeta function, a pole of T(sJ2) at s « 0 with 
residue (+2) and a pole at s - 1 with residue (-2) due to T(i-j/2), and find 

1 < c (A9) 
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To obtain this expression, note that the contribution from the pole at s = 1 cancels 
with the term in IJD in eqn (A8). We have also used the Riemann relation 97 
linking £(s) with ((1 — s), a corresponding relation 96 for fi(s) and the duplication 
formula for the gamma function. 97 If we now use eqn (A6) again the last integral 
can be easily re-expressed as a sum, and since ((i) = -1.460, fid) = 0.6677, 96 we 
have 

2n2>e 2 r 0.6206^ • 1 f 2nD .— 2 — 5 H iL 

* - -srL 1 — oh>£. pro? ■» r^ n + " 11 <A10) 

Note that if 2> ^ ^/a, the sum converts with extreme rapidity and can be 
ignored, so that to sufficient approximation 

"ei * -HI ~l (All) 

For small D/yfa 9 the corresponding expansion follows directly from eqn (A8) and 
gives 

where C(i) = 2.612, JB@) = 0.864. 

These techniques can also be developed further to deal with Stern layer problems. 
Note that for D ~ yfa, the " effective M capacitance distance D is reduced by a 
factor (1 -0.62yfa/D). If image effects are included, it can be shown that (when the 
zwitterions are immersed in a high dielectric medium a 80, adjacent to a hydrocarbon 
surface e = 2), image effects can virtually double the electrostatic energy so that the 
two effects may partially cancel out. Nonetheless the observation that a 0l is reduced 
from the intuitive capacitance form is of some interest, as it is known * that for ionic 
micelles use of Debye-Huckel theory (at low salt equivalent to the capacitance form) 
gives results for the repulsive free energy too large by a factor of 2. 



APPENDIX B 

(1) MINIMUM FREE ENERGY FOR VESICLES 

We outline here the steps necessary to derive the results quoted in section 7. 
(a) spherical vesicles: no packing constraints 

Consider first a spherical vesicle comprising N amphiphiles. The free energy of 
the vesicle is 

where R u ^il Rt> n 2 refer to outer and inner surfaces. When packing is ignored, the 
only constraints are those due to conservation of volume and number, viz 

Nv - y(H?-*S) (B2 

N=n x +n 2 . (B3) 
Since N is fixed, there are only two variables which can be taken as n u R x . Then 
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n 2 and R 2 are determined. We first minimise G with respect to n u taking R u R 2 
constant. This gives 

dG 

drit 



/%^Y - 2ne2D ^ 2RiR ^ 



(B7) 



This equation can be solved by iteration as follows : clearly, for R lt i? 2 > D there is 
only one solution. The first approximation can be found by dropping terms in D 
inside the square brackets. This gives 



or a t = a 2 - a 0 . To find the second approximation to the solution of eqn (B7), 
we find solutions R u R 2 of the equations 



Substitution of these values into eqn (B7) then gives the average area a to second 
approximation 



(B12) 



3 



U e +Z>) K 2 (* 2 -I>)j V ' 

Physically this condition is equivalent to the statement that there is no potential 
difference across the vesicle in a situation of minimum free energy. Using eqn (B3) 
and (B4) we have 

We now minimise with respect to the remaining variable R u using dR 2 /8R x = 
R\IR h which follows from eqn (B2). The condition for a minimum is 

^<m<G) = 0 (BQ J 

or, after a little algebra 




4jr(K?+R|) = jV flo ; i 



^(K?-^)-^. (B9) I 



These approximate values can then be substituted inside the square brackets of 

eqn (B7) to yield a better value for Rl+Rl on the left side. Writing t = R t ~R 2 , ?' ; 's| 

the solution of eqn (B9) is 



■a 
J 

f 
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Optimal areas a u a 2 are now given as 



AnR\ f D(R 1 + 3R 2 ) } 
ej= _ =B „| 1+ _L_^. (8,3, 



The minimum free energy is now from eqn (B5) 



- 2tfya 0 -4rtyZ>*, (B14) 
where / is given by eqn (BIO). 

(6) CRITICAL RADIUS 

Let / c be the maximum extended length of a lipid molecule. As N decreases we 
eventually reach a critical radius below which packing constraints become important. 
This radius R\ is determined by 

^ = 4 ^ = $C*?-(^-0 3 ] (B15) 
with a t given by (B13). To a sufficient approximation a x « a 0 (B15) becomes 

*c = -^Tl + -^5 +f ( -TS Y+ • • -1 (B16) 
Eqn (B15) determines R c x and is a special case of eqn (5.4). 

(c) PACKING CONSTRAINTS 

Now suppose that N < N e . We wish to find the configuration of such a vesicle 
with minimum free energy. The additional constraint is now eqn (B15), and for 
fixed N> 1 € the free energy involves only a single variable. Note now that we must 
have a potential across the membrane. It is convenient here to rewrite eqn (Bl) in 
the form 

G = «i^i+^)+^2+0 (B17) 

where we have ignored curvature corrections to the electrostatic contribution to a 
first approximation. A complete solution to the problem of minimising this function 
is very difficult and we proceed as follows: we first show that da 2 jda x > 1. This 
follows from the constraints 

y[^-(*i-U 3 ] = «,»^^-= —(2^-0 

dn 2 4nl c ,„_ 
„ I+ „ 2 = N => -f = — — K2Ri-Q 
oK x v 
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- c (21?i-/c)| 



= a, => 



L 2 




Hence 



n 2 2 dR x \n 2 Ri n \» J 
Jj^^^lR, -/I 

d^-SRjdR, *IR\ L^(2J? r !i 

1 «i» J 



Now from eqn (B15) 

-^ (2J?1 - U= "»7"V 2Kj 2J?J V 2RJ 

Hence (B18) becomes 

£22l > ^/^A (B20) 
Tvmcafly for egg PC, * a - 100 A, * 2 ~ 70 A, h ~ 17 A, Ife/fel > 30- Con- 

therefore in G. Thus we can conclude that wnen ir is neur i» ^ 
fl2 ~ fl 0 , and after rewriting (B17) in the form 

G = 2Nyao+ «J^^^"j2^^ 2 ' (B2D 
we may ignore the term in n 2 . Thus G = 2^o+Wi(l-«o/*i)*. From (B15), 

jTyRj - 1 ~rX R\) 
G«2Nya 0 +-^-^l-^J 

The corrections to this result due to curvature can, to the lowest approximation, be 
computed ignoring packing, and added to (B22). 

(2) TUBULAR VESICLES 

Consider now a tubular vesicle. We shall show that the free energy per lipid of 
snch^ratfemotyts^ater than that for a bilayer, so that tubular vesicles w.11 not 



(B16) 
Hence 
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form, being disfavoured both energetically and entropically. To do this, it is not 
necessary to deal with the more difficult problem which includes complications due to 
packing. In the capacitance model for the repulsive forces, the free energy per unit 
length of a tubular vesicle is 

where R u R 2 are outer and inner radii, n u n z denote number of head groups per unit 
area. For a given aggregation number N per unit length, the constraint equations 
corresponding to eqn (B2), (B3) are 

N m n x +n 2 > *{R\-Rl) = Nv. (B24) 

Again we take n u R x as the two independent variables, and first minimise with respect 
to n x and fixed R x remembering that n x +n 2 = N. This gives 

mb,,g = 2^ 1+Sl ) + ^h(^±i). (B25) 

Now, minimising with respect to J?, we find 

| min = min* t |min ni ^jl = 2jiy(* 1 +i? 2 )+ 



e 



(B26) 



It is now straightforward to demonstrate that the optimal free energy per lipid GfNL 
(min) satisfies 



G 



W>{«+^1+j^]>2^ (B27) 

If the packing constraints are included the free energy per lipid would be even larger 
than 2ya 0 , the optimal energy for a bilayer. Hence we conclude that tubular vesicles 
can not exist, Le. 9 that planar or spherical bilayers are the preferred assembly. 
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